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^ \ Abstract: We construct a five dimensional supersymmetric SO(10)xD3 grand unified 

model with an S 1 / (Z 2 x Z' 2 ) orbifold as the extra dimension. The orbifold breaks half of 
D ! the supersymmetry and breaks the SO(10) gauge symmetry down to SU(4) C x SU(2) L x 

SU(2) R . The Higgs mechanism is used to break the remaining gauge symmetry the rest 
of the way to the Standard Model. We place matter fields variously in the bulk and 
on the orbifold fixed points and the resulting massless fields are mixtures between 
these brane and bulk fields. A chiral adjoint field in the bulk gets a U(l)x vacuum 
expectation value, resulting in an X-dependent localization of the bulk matter fields and 
the Standard Model Higgs field. This Higgs field localization allows us to simultaneously 
explain the hierarchies m u < rrid and m t ^> m&. The model uses 11 parameters to fit 
the 13 independent low energy observables of the quark and charged lepton Yukawa 
matrices. The model predicts the values of two quark mass combinations, ^ and 
mdirismb, each of which are predicted to be approximately la above their experimental 
values. The remaining observables are successfully fit at the 5% level. We note that 
this 5D theory, as formulated, has problems retaining the Pati-Salam Yukawa symmetry 
relations. A simple 6D fix, which preserves the 5D results, solves this problem. 
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1. Introduction 

The Standard Model (SM) of particle physics has had spectacular success in describ- 
ing the strong, weak, and electromagnetic forces of nature in terms of gauge theory. 
It does, however, leave several issues unexplained. Among these issues: First, in the 
Standard Model the three gauge coupling values are arbitrary and unrelated, but ap- 
pear to almost unify at a high scale. The assumption of supersymmetry improves this 
unification. Second, the fermion charge assignments under the gauge symmetries are 
arbitrary. Third, recent experiments have shown that neutrinos have masses, but they 
are massless in the Standard Model. Fourth, in a universe initially balanced between 
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matter and anti-matter, interactions which violate baryon number are required to ex- 
plain the apparent baryon asymmetry observed in the universe. Baryon number is 
conserved by the renormalizable terms of the Standard Model. Fifth, the weak and 
Planck scales are 16 orders of magnitude apart, causing a naturalness problem in the 
mass of the Higgs boson. While this is not an exhaustive list of problems with the 
Standard Model, all of these issues are addressed by four dimensional supersymmetric 
grand unified theories [4D SUSY GUTS] 0, g § |,|, §. 

In a SUSY GUT, the three gauge couplings meet near the GUT scale Mq ~ 
3 x 10 16 GeV, and above this scale an enlarged symmetry group with a single gauge 
coupling governs physics. Within this paper, we will take SO(10) as our gauge group. 
Remarkably, the charges of the SM fermions are exactly those which come from a 16 
of SO(10), including the arbitrary (in the SM) U(l)y hypercharge quantum numbers. 
Moreover, a 16 dimensional representation requires the addition of a sterile neutrino 
which can lead to tiny neutrino masses through a See-Saw mechanism. B — L is a 
gauge symmetry within SO (10) which is broken at a high scale in the unified theory, 
thus making Baryogenesis possible. Finally, the presence of supersymmetry near the 
weak scale gives a natural explanation for a Higgs mass at that scale. 

As experiments become more accurate, we are able to test whether GUT theories 
work in their minimal forms. There are several indications that the most simple GUT 
theories cannot work. 

• With the definition of the GUT scale M G : ai(M G ) = a 2 (M G ) = «gut, the as- 
sumption of exact unification (a 3 (M G ) = «gut) leads to a prediction for a 3 (M z ) 
which is larger than the measurements. This indicates that a negative threshold 
correction to a 3 must be present at the GUT scale: £3 = a3 ( MG )~ aGUT ~ —0.04 
0. 

• Since the fermions are combined into fewer multiplets of the unifying group, 
GUT theories in general give relationships between the fermion masses at the 
GUT scale. In SO(10), third family unification can be accomodated (m T (Mc) = 
mbiMc) = mt(Mc)), but the first and second family masses do not unify at Mq- 
The Georgi-Jarlskog relation: m s /m^ = |m;,/m r ||, which had been known 
to work very well, holds less well as measurements of the strange quark mass 
decrease. 

• GUT models place the weak Higgs doublet (s) into larger representations which 
include new color triplet Higgs fields. It is difficult to make the color triplets 
heavy enough (near GUT scale mass) to avoid rapid proton decay while keeping 
the weak doublets light (weak scale mass). The doublet-triplet (DT) splitting 
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problem is related to the method by which the GUT symmetry is broken to the 
Standard Model, and the models which give proper DT splitting can be quite 
complicated. 



Orbifold GUTS have been considered extensively for the last five years pL 10, 11 



|12j , jig , |I4[ , |15| , |16| , [H], [L8], [19] and provide answers to some of the above issues. If the 
compactification radius is slightly smaller than the cutoff scale, the heavy Kaluza-Klein 
(KK) modes can give compactification scale threshold corrections necessary for exact 
gauge unification at the cutoff scale |p0| . If the Standard Model Higgs is a bulk field, 



and if the triplet Higgs has twisted boudary conditions under the orbifold it will acquire 
a mass at the compactification scale, thus solving the DT splitting problem naturally. 
The same mechanism breaks the GUT gauge symmetry and gives mass to gauge bosons 
outside of the Standard Model. 

5D SU(5) models have a SM brane on which all of the interesting GUT mass re- 
lations are absent. The first and second family matter fields usually feel this brane 
and so there are no relations between quarks and leptons in these two families. This is 
unattractive since the models lose predictivity. For this reason, we choose to concen- 
trate on SO(10) models which give more possible avenues of symmetry breaking to the 
Standard Model and allow for GUT mass relations. 

5D SO(10) SUSY GUTS were first considered in Dermisek and Mafi [[L8[], and a 



setup which gives gauge coupling unification was constructed by Kim and Raby 
The extra dimension is an S l /(Z 2 x Z' 2 ) orbifold, a line segment with endpoints which 
are fixedpoints of the orbifold Z 2 symmetries. 5D N=l SUSY is broken to 4D N=l 
by the first Z 2 , while the second Z 2 breaks SO(10) to SO(6) x SO(4) (= SU(4) x 
SU(2) x SU(2) = Pati-Salam gauge group). Further breaking to the Standard Model 
is provided by the Higgs mechanism on one of the fixed points. One fixed point, the 
SO(10) brane, is invariant under the first Z 2 and has SO(10) gauge symmetry. The 
other fixed point, the Pati-Salam (PS) brane, is invariant under the second Z 2 and only 
has Pati-Salam gauge symmetry. The bulk has the full SO(10) gauge symmetry. Gauge 
coupling unification works if the compactification scale is approximately 10 14 GeV, the 
cutoff scale is approximately 10 17 GeV, and if the Higgs multiplet lives in the bulk. A 
virtue of the model is that SO(10) or Pati-Salam exists after the orbifold breaking and 
we retain some of the Yukawa GUT relations. 

We note that there is a problem with the five dimensional formulation of our theory. 
In order for the gauge couplings to unify at the cutoff scale, the PS brane fields breaking 
Pati-Salam down to the Standard Model must have cutoff-scale VEVs. Therefore, in 
general one would expect all aspects of PS symmetry to be broken at this high scale, 
and that the lower-energy theory should not be expected to exhibit PS symmetry. 
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Nevertheless, we wish to retain some of the PS symmetry below the cutoff scale in 
order to preserve the PS Yukawa coupling relations. This is possible if there is an 
additional extra dimension in which the PS brane is separated from the PS-breaking 
VEV fields. This 6th direction will solve this problem if its length scale is only slightly 
larger than the cutoff scale. The effective 5D theory below the energy scale of this 6th 
direction will not be affected in any other way. For this reason we choose to perform 
our analysis in the 5D theory, and will assume that PS Yukawa relations are still valid 
below the cutoff-scale. We discuss the 6D fix to the PS-breaking problem in Section |5|. 

Attempts to explain the mass hierarchy between the different families of fermions 
have often centered on flavor symmetries. There are models based on abelian horizontal 
U(l) symmetries where the smallness of certain couplings is explained by the suppres- 
sion given by high powers of U(l) breaking fields. The explanations can at most be 
qualitative as there are order one coefficients which are undetermined from the U(l) 
symmetry. We therefore choose to concentrate on non-abelian flavor symmetry. With 
three families, the largest possible symmetry would be SU(3). However, the order one 
top Yukawa coupling would badly break the SU(3) symmetry, and so instead we choose 
to concentrate on an SU(2) symmetry between the first two families. This symmetry 
can explain the absence of flavor changing neutral currents in supersymmetric theories 
and can relate unknown order one coefficients between different families. In string the- 
ory, global symmetries are thought to be generally broken by quantum gravity effects 



21j |22| , |23| , and so our flavor symmetry should be a gauge symmetry. However, the 



breaking of a continous gauge symmetry like SU(2) gives unwanted flavor changing 
neutral currents from the D-term contributions [|4], |25| . For this reason, we assume 
a flavor symmetry relating the first and second families based on D3, a discrete sub- 
group of SU(2). As a discrete gauge symmetry, D 3 has all the virtues of SU(2) and can 
avoid the problem related to continuous gauge symmetry |2"B| . More information and 
references on family symmetries in grand unified theories can be found in the reviews 



In SUSY SO (10) models with flavor symmetries, the most difficult thing to explain 
is why m u < m d while m t m b . In SO(10) models, we usually have m t /mi, = 
Y t /Y b tan/? ~ tan/3, and the heaviness of the top quark is explained by a large value 
for tan/? ~ 50. In these models, it is natural that Y u /Yd ~ Y t /Yb, which leads to 
TTiu/nT'd ~ 50. To fit the data, an unusually small Yukawa coupling is needed for the up 
quark compared to the down quark. This is difficult to implement in SO(10) models. 
Extra-dimensional theories provide a nice tool to suppress or enhance fermion masses, 
as the size of a mass can be influenced by the localization of the matter and Higgs fields 
within the extra dimension. Others, notably Arkani-Hamed and Schmaltz in ||30|| , have 
used localization in an extra dimension to explain the fermion mass hierarchy. Our 
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model, on the other hand, retains some of the attractive Yukawa relations present in 



4D GUT models which are not present in [30 



We explain the ratios m u /m^ < 1 and m t /m b 3> 1 simultaneously by the quasi- 
localization of the Higgs field in the bulk by a kink mass. 1 We give a vacuum expectation 
value (VEV) in the U(l)x direction to a scalar adjoint field £ contained in the 5D 
vector multiplet. S is odd under both Z 2 parities, and we choose the VEV so that 
its contribution to the D-term VEV is only at the two branes. Supersymmetry can be 
preserved if there are fields on each brane which get VEVs in the right-handed neutrino 
direction to counter the contribution of (S) to the D-term VEV. The VEV of £ acts 
to give a U( Independent mass to the hypermultiplet fields in the bulk, in particular 
the Higgs field. Because H u and H d carry opposite U(l)x charges, they are localized 
towards different branes. Thus if the 1st and 3rd families get their Yukawa couplings 
on opposite branes, we can naturally have m u /m^ < 1 and m t /m b > 1. 

The easiest way to ensure that the 1st and 3rd families get their Yukawa couplings 
on opposite branes is to restrict them to these branes. By proton decay constraints, 
the 1st family cannot reside on the SO(10) brane. Thus we choose the 1st (and by 
family symmetry, the 2nd) family to be on the PS brane, and the 3rd family to be on 
the SO(10) brane. We then choose the sign of the kink mass VEV to be such that 
H u is localized towards the SO(10) brane and H d is localized towards the PS brane. 2 
Communication between the first two families and the third is provided by mixing with 
bulk fields. After requiring that the 1st and 3rd families be brane fields, the location of 
the Higgs and all matter fields is entirely determined from gauge coupling unification, 
proton decay constraints, D 3 family symmetry between the first two families, and the 
action of the kink mass to explain 1st and 3rd family mass ratios. 

To summarize: In this paper, we construct a 5D SO(10) supersymmetric orbifold 
GUT model in which the 1st and 2nd families reside on the PS brane and transform as 
a doublet under the D3 family symmetry, the 3rd family is located on the SO(10) brane, 
and a kink mass localizes the MSSM Higgs doublets to opposite sides of the bulk. The 
model explains the 13 independent quark and charged lepton masses and mixing angles 
in terms of 11 parameters. The two predictions are — (Mz) = 0.0037 ± 0.0006 and 



lr rhe extreme limit is discussed in where H u is on the brane and Hd is in the bulk such that 
rat 3> nib for order one tan (3. 

2 These arguments do not rule out the case of all 3 families in the bulk, or some fields on branes 
with others in the bulk. However, if the 3rd family is a brane field, then the 2nd family must also be 
a brane field or the ratio ^f- is made too small by volume suppression to fit the data. It is possible 
to decrease the suppression from the volume factor by either placing only half of the 2nd family in 
the bulk and/or by altering the GUT threshold correction £3. We made some attempts with these 
approaches but were unable to produce viable theories. 
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mdm s mb(Mz) = (10.7 ± 5.0) x 10 5 MeV 3 , both of which are approximately la above 
the experimental values. 

We lay out the remainder of this paper as follows. The data is summarized in 
Section ||. The basic setup is explained in Section [3[ This includes background material 
and our choice of superpotential and Yukawa matrices. Section f| contains the analysis 
of our data. This part is broken into two subsections: an analytic section (Section |4.1|) 
giving predictions and approximate relationships between the model parameters and the 
data, and a numerical section (Section f4.2| ) with more precise results. The possibilities 
and advantages of placing our model in a 6D framework are given in Section |5]. Finally, 
we present a summary and discussion of our model in Section ^| We also provide 
necessary information on D3 family symmetry in Appendix [A| and a determination of 
massless state wavefunctions in Appendix |B|. Finally, in Appendix we consider a 4D 
SO (10) SUSY GUT model |32j lifted up to 5D. In this model the smallness of the up 
quark mass is explained by an approximate left-right (LR) symmetry rather than by a 
kink VEV. 



2. Data 



In this section we tabulate the low energy data used in our analysis. First the data as- 
sociated with the Cabbibo-Kobayashi-Maskawa (CKM) matrix. These CKM elements 



we take from 33 



IK, 
\v cb 



v. 



cb 



\Vtd\ and sin 2/3 we take from []3~4 |. 

\Vtd\ 
sin 2f3 



0.2240 ±0.0036 
(41.5 ±0.8) x 10" 3 

0.086 ±0.008 



0.0082 ±0.0008 
0.739 ± 0.048 



We use |35[ for J and Ek- 

J = (3.0 ±0.3) x 10~ 5 
e K = (2.282 ±0.017) x 10~ 3 

Next the low energy quark mass observables. We use 



Q = 



1 - 



1 »d 



(2.1) 



(2.2) 



(2.3) 



(2.4) 
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and rris/md from |3B|] with doubled errors on m s /m^, the unquenched lattice QCD 
result with nj = 2 for m s [^7] with doubled errors, and we take the quenched lattice 



QCD result for m c with 10% error from |33| and double the error. We use the bottom 
quark mass from and the top quark pole mass from the CDF and D0 Collaboration 



Q 

m s 
m d 

m s (2GeV) 
m c {m c ) 
m b (m b ) 
M t (pole) 
m t {m t ) 



22.7 ±0.8 
18.9 ±0.8 x 2 

89 ± 11 x 2MeV 
1.30 ±0.15 x 2GeV 
4.22 ±0.09 GeV. 
178.0 ± 4.3 GeV 
169 ± 4 GeV 



(2.5) 



In our analysis, we will calculate our observables at the energy scale Mz, except for 
the top mass for which we will use m t (m t ). Therefore we need to convert our data to 
values at Mz- Define the running parameters: 



At two loops in QCD we find 



m[Mz} 

mi(mi) 

mj(Mz) 
mi (2 GeV) 



Vc 

rju 



0.56, 



for % 



for i 



Vb 



= c, b 
= u, d, s. 

0.69, 
0.65. 



(2.6) 



(2.7) 



The lepton masses at Mz from |4^| are 



m e 



0.48684727 ± 0.00000014 MeV 
102.75138 ±0.00033 MeV 
1746.7 ±0.3 MeV. 



{21 



Also from i40fl , we take the gauge couplings at Mz'- 



a x {M z ) 
a 2 {M z ) 
a 3 (M z ) 



0.016829 ±0.000017 
0.033493 ± 0.000059 
0.118 ±0.003 



(2.9) 
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We use the mass of the Z-boson to set the weak energy scale: ME 

M z = 91.1876 ± 0.0021 GeV 



(2.10) 



In our model, we do not calculate the effects of electro-weak symmetry breaking. We 
assume that it happens properly and gives the correct weak scale masses. For this 
reason, we use a calculated Higgs VEV at M z . 



7i^la 1 {M z ) + a 2 {M z ) 
246.41 ±0.17 GeV 



(2.11) 



Additional information is needed to calculate Ek in our model. We use the CKM 
parametrization- independent formula 



a 



£ B r 



I | = 



Lin [ S c a c + Im ( S t a t + 21m J S,,cr rl 



a. 



where 



□i = v ud v: s v is v* d 



a. 



(2.12) 



(2.13) 



and the various Si functions are Inami-Lim functions f4l| . The Oi are 0(1) factors 
which we take from Battaglia [33| 



1.32 a t = 0.57 a 



ct 



0.47. 



(2.14) 



C e is a ratio of well known low-energy observables. From Battaglia j33 , this combina- 
tion is 



a 



3.837 x 10 4 . 



(2.15) 



We neglect the error associated with C e since it is small and will not affect the theo- 
retical error in A large part of this theoretical error comes from Bk for which we 
use 



33 



B K = 0.86 ±0.06 ±0.14 
~ 0.86 ±0.15 



(2.16) 



where in the last line we add the statistical and systematic errors in quadrature. 
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3. Setup 



3.1 Background 

We consider a five dimensional supersymmetric SO(10) GUT compactified on an S 1 / (Z 2 x 
Z' 2 ) orbifold where S 1 is described by y G [0, 27rR]. The first orbifold, Z 2 , under which 
y — > —y, breaks 5D N=l supersymmetry (4D N=2) to 4D N=l. The other orb- 
ifold, Z 2 , under which y — > — y + ttR, breaks SO (10) down to the PS gauge group 
SU^)^ x SU(2) L x SU(2) R . The fundamental domain of the y direction is the line 
segment y G [0, ttR/2]. SO(10) gauge symmetry is present everywhere except at the 
point y = nR/2, which only has Pati-Salam gauge symmetry. Hence, we call the two 
inequivalent fixed points the "SO(10)" (y = 0) and "Pati-Salam" branes (y = nR/2) 
where each fixed point is a three-brane (3+1 dimensional spacetime). The Higgs mech- 
anism on the PS brane completes the breaking of the PS gauge symmetry to the SM 
gauge group. 

The fields which live in the five dimensional space between the branes (known as 
the "bulk" ) are even or odd under the orbifold parities. We denote the parity of a given 
field under these orbifolds by two ± subscripts on the field with the first corresponding 
to Z 2 and the second corresponding to Z' 2 . This part of the setup, including the orbifold 
structure, field parities, and supersymmetry and gauge symmetry breaking are based 

For more details, please see these references. 



on work done in 18 and 20 



We wish to relate some of our fields by a family symmetry, D3, under which the first 
and second family fields form a doublet. Other fields within our model will be in various 
representations of D3 which will affect the structure of the Yukawa matrices we generate. 
We take the family symmetry to be independent of the orbifold symmetries. A brief 
summary of the D 3 group is provided in Appendix [A], where we give the information 
necessary to understand the family representations and couplings used in our model. 
For more information, see the Appendix presented in [ |32|j . 

The 5D supersymmetric vector multiplet V = (Am, Ai, A2, S) contains a 4D vector 
multiplet V = (A^, Ai) and a 4D chiral adjoint <\> = ((£ + iA 5 )/^/2, X 2 ). For a generic 
hypermultiplet H = (h,h,ip,ip) which breaks up into the 4D chiral multiplets $ = 



(h,ijj) and $ = (h,ip), we use Arkani-Hamed et. al. [42| for the 5D action: 



S D / d A xdy 



d A e 



(3.1) 



d 2 9 $(m + d y 1=4>)^ + h.c. 

\/2 



The mass m refers to the mass of the 5D fields, which we set to zero for all hypermul- 
tiplets. We take the field E (c </>) to get a VEV in the X direction of SO(10). Because 
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of the coupling in the action: 



S D d 4 xdy d 2 9 



7T 



(3.2) 



VEV 



this VEV will generate V-dependent masses mx for the bulk chiral fields coming from 
hypermultiplets in the theory. The XJ(l)x of SO(10) is contained inside Pati-Salam, and 

the chiral adjoint for PS has ( ) parity. Therefore the mass mx we have generated 

also has this odd-odd parity under the orbifold. For this reason we call it a kink mass. 
The D-term for this theory is 



D = -(d y X+i[Z,A 5 )), 



(3.3) 



hence a VEV of E is potentially dangerous since it could create a D-term VEV. To 
avoid this, we take (E) to be flat in the bulk with discontinuities at the branes to be 

consistent with its ( ) parity. These discontinuities generate D-term VEVs on the 

branes, but we can choose to add brane fields (16kmk on the SO(10) brane, Xkink an d 
X c kink on the PS brane) with VEVs in the v c and v c directions to cancel these effects 
from <9 y (E), leaving us with a D-flat theory. 3 More concretely, if the sterile sneutrino 
component of 16kink Kink) an d the charge conjugate sterile sneutrino component of 
X c kink (^ c kink) get their VEVs, the D-flat condition for the U{l)x subgroup is given by 



5|Kink)l^(?/)-5K^kink)| 



(3.4) 



As ^ ink and z/ c kink carry only U{l)x charges, they do not appear in the D-term of the 
standard model gauge group. The D-flat condition is satisfied for 



(E) = E e__(j/), 



I Kink) P 



IKkink)| 2 



2Ec 



(3.5) 
(3.6) 



E- + and e. 
analysis. 



are step functions on the orbifold. Both of these will be used later in our 



■(v) = 



-1 for ye [-ttR,0] 
-1 for ye [0, ttR) 

' +1 for y E [-irR,-!f] 



-1 for y E [-^,0] 
-1 for y G [0, 
-1 for y E [^f^R] 



ttRi 
2 J 



(3.7) 



(3. 



3 This technique is described in W2| and used in 
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We choose to parametrize the kink mass as follows: 

m x (y) = e_(y)X^ = *-(v)*^ (3.9) 
C ~ 2 ' 

( will turn out to be a useful dimensionless parameter in later analysis. 

If our £/(l)x-breaking fields were to get VEVs as large as the cutoff scale, we could 
not keep SO (10) or Pati-Salam as our symmetries on the branes. However, it turns out 
that £ ~ 2 is needed to fit the observed physical quantities, and so the corresponding 
VEVs are 



Kink) 



5nRg 2 




Therefore, all the VEVs necessary to give a kink mass are around the compactification 
scale and the breaking effects are suppressed at least by M c jM* ~ 10~ 2 or 1CT 3 P0 



For the most part, the U(l)x breaking effects come from the kink profiles of bulk fields 
which are calculable and are exponentially proportional to the U(l)x charges. It is 
a novel example of obtaining large (order one) symmetry-breaking with a very small 
symmetry-breaking order parameter S. We stress that the VEVs of E, i/£ ink and ^ c kink 
do not spoil the symmetries on the branes. 

On the other hand, on the Pati-Salam brane gauge coupling unification requires x° 
and x c t° get VEVs of order the cutoff scale (x c and \ c are different from x c kink an d 
Xkmk)- This is a serious problem for our theory. There are order one corrections in the 
Kahler potential and Pati-Salam symmetry is badly broken in the canonical basis: 

K= + (3.10) 

where C is of order one, (x c ) = (x c ^) ~ M*, and where if> is a generic Pati-Salam brane 
field. This problem can be solved in a geometric way if there is a sixth dimension along 
which x c an d ip are separated. If this sixth dimension has a length scale R2 slightly 
larger than the cutoff length scale, we can simultaneously retain the desired PS-brane 
Yukawa relations and gauge coupling unification. The 6D setup will be discussed in 
detail in section |5|. Until then, we choose to concentrate on the effective 5D theory 
below this 6D scale. The 5D analysis of the Yukawa matrices will not be quantitatively 
affected by the addition of this extra dimension. 

We now briefly summarize the gauge unification results given by Kim and Raby 



20[ . A 5D gauge theory is nonrenormalizable and gets large corrections at the cutoff. 
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Corrections to gauge couplings, however, will be the same for all couplings unified into 
the larger gauge group. These corrections will affect the absolute values of the gauge 
couplings, but not the differences. Further, if the gauge symmetry is broken only by 
orbifolding or by Higgs mechanism on the branes, the differences in the couplings will 
have a logarithmic, calculable running. 

The states which affect the differential running are the bulk vector multiplet V 
and the bulk Higgs hypermultiplet 7i. 4 The placement and number of complete matter 
multiplets does not affect gauge coupling unification, since matter multiplets (16s of 
SO (10)) act equally across the three gauge couplings and cannot affect the coupling 
differences. Those states in the theory outside of the MSSM have twisted orbifold 
boundary conditions and so have masses at the compactification scale (M c ). For en- 
ergies below M c , the theory is the MSSM. The effects of running between M c and 
M* (the cutoff scale), including the Kaluza-Klein (KK) towers, are taken as threshold 
corrections at M c . Without these threshold corrections, it is known that the MSSM 
unifies around Mq ~ 3 x 10 16 GeV with a coupling of «gut ~ 1/24 and a GUT-scale 
threshold correction for a 3 of £3 ~ —0.04. Assuming unification in the orbifold theory 
at the cutoff scale M* and that the PS breaking Higgs VEV is of order the cutoff scale, 
we can solve for M* and M c in terms of the 4D GUT parameters Mq, ogut, and £3. 
This leads to M* ~ 10 17 GeV and M c ~ 10 14 GeV f20|. 

The matter field locations are constrained by proton decay: 



Matter fields on the SO(10) brane 

There are gauge bosons within SO(10) which mediate baryon (B) and lepton (L) 
number-violating interactions. All of these are outside of the Standard Model 
and hence have masses of order M c or higher. After integrating out these fields 
(and their KK modes), we get dimension six operators which violate B and L for 
any matter multiplets on the SO(10) brane. These operators are suppressed by 
1/M^. Given M c ~ 10 14 GeV, current bounds on proton decay rule out models 
which have these operators for the first and second families. Thus only the third 



family can reside on the SO(10) brane [fill. Dimension 5 proton decay operators 



4 Whcn considering differential running of the gauge couplings, a Higgs hypermultiplet in the bulk is 
effectively the same as a 4D 10 of SO(10) with light Higgs MSSM doublets and heavy Higgs triplets of 
mass M c . This setup admits gauge coupling unification as shown by Kim and Raby |2(J. In particular, 
see the calculations leading to equation (3.13) of that paper. Effects from brane Higgs doublets would 
be felt up to Af» and would tend to inhibit unification since they drive the couplings apart rather than 
together. 

5 With dimension six operators for the third family, mixing between this family and the first two can 
induce proton decay. Assuming that the mixing is of order \ V u b\ or | V c b\ and that the gauge bosons have 
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vanish since the color triplet Higgs states obtain off-diagonal mass with triplets 
in 10 and these states do not couple to matter. 

Matter fields on the PS brane or in the bulk 

Pati-Salam gauge symmetry does not relate the left-handed fields if) ((4,2, 1) in 
PS) to the right-handed fields ip c ((4,1,2) in PS), and we do not get baryon 
number-violating dimension six operators after integrating out the heavy gauge 
bosons. Therefore, any matter fields can be on the PS brane as long as the PS 
breaking scale is not extremely low. This scale in our theory is M* ~ 10 17 GeV, 
and proton decay is not a problem here. 

In principle, we can consider higher dimensional operators with derivative inter- 
actions c?5 = d/dy. Because the coefficients of the higher dimensional operators 
are not determined from the theory we cannot calculate the proton decay rate 
from these operators accurately. However, we can get a bound that is consistent 
with our setup by assuming unknown coefficients to be order one. See Kim and 



Raby for more details 20 . 



Proton decay constrains the first and second families to reside either in the bulk or 
on the PS brane, but does not constrain the location of the third family. We choose to 
place the third family on the SO(10) brane. Given this placement, the second family is 
forced to reside on the PS brane. If it were in the bulk, volume suppression would give 
m /i/ m T ~ 10 -3 , which is far too small to fit the data. Our D 3 flavor symmetry places 
the first and second families together into a doublet, so we place both on the PS brane. 

Let us summarize the basic setup. 

• Gauge symmetry : SO(10) in the bulk and at y = 0, Pati-Salam at y — ttR/2. 

• Higgs fields come from 10 dimensional hypermultiplets in the bulk. 

• 3rd family matter fields are on the SO(10) brane. 

• 1st and 2nd family matter fields, a doublet under D 3 , are on the PS brane. 

• Kink mass localizes the hypermultiplets through their X-dependence. 

mass at M c , naive calculations using formulae in ftifl put the proton lifetime many orders of magnitude 
above current limits, since this leads to an effective gauge boson mass of order M c /(V c bV u b) ~ 6 x 10 17 
GeV. 
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Table 1: Bulk fields 



Field 



PS Symm D 3 Symm 



10 



10 



16 



16 




H 

TT C . 



+ 



) 




3.2 Yukawa Matrix 

This section introduces the fields and superpotential of our model. Here we calculate 
the Yukawa matrices associated with the massless fields corresponding to the Standard 
Model fermions. 

First we introduce the matter content of the theory. In the bulk, we have two 
5D N=l hypermultiplets which transform as 16s under SO(10). These fields form a 
doublet (2a) under the D3 family symmetry. In our Lagrangian, we show the D3 doublet 
structure of these fields by a subscript a, with a = 1, 2. For information on how doublets 
and other D 3 objects couple, please see Appendix [X[ We also have a hypermultiplet 10 
of SO(10) which is a singlet (1a) under the family symmetry. These fields are listed 
in Table |I| with their parities, PS gauge symmetries, and family symmetries. On the 
SO(10) brane we place a single 16 of SO(10), invariant under D 3 , and an SO(10) gauge 
singlet which is a doublet under the family symmetry. These fields are listed in Table 
0. The PS brane has a more complicated set of fields. There are the fields ip and ip c , 
each doublets under D3, which transform as the two halves of an SO(10) 16. N c and iV c 
are also family doublets and carry charge under PS in order to allow mixing between 
ip c and N c . There are various other fields 0, A, A 3R , u c , and u c which will get vacuum 
expectation values. These fields are listed in Table |3| with their PS charges and D3 
family symmetries. 

We choose to break the superpotential into two pieces: W = W\ + W 2 - The first 
contains terms leading to interactions between the matter and Higgs field (H = H ++ ): 



W, 



-A3I63IOI63 <%) 
+ \\ a ip a HipaA + Xfsij a HN c J a + 2\^ a HA m (C+)J 5(y 



irR 



2 



(3.11) 
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Table 2: SO(10) Brane fields 
Field PS Symm D 3 Symm 

(1,1,1) 2 A 

Table 3: PS Brane fields 



Field PS Symm D 3 Symm 





(4,2,1) 


2 A 




(4,1,2) 


2a 


N c 


(4,1,2) 


2a 




(4,1,2) 


2 a 


4> 


(1,1,1) 


2a 


A 


(1,1,1) 


Ib 




(1,1,3) 


1a 




(4,1,2) 


1a 


~ c 


(4,1,2) 


1a 



W 2 contains mass terms: 

W 2 = TE a (d y -m x )16a + 10(d y -mx)10 (3.12) 
+ [2a0 a T6 a 16 3 ] 5{y) 

_ 7]-/? 

+ [2 V (^ ++ ) a Xa + N° a (aft + bo^K)} 8(y - — ) 

Sufficient factors of the cutoff scale M* should be placed in the superpotential terms 
so that the couplings A 3 , X a , \p, A 7 , a, r], a, and b are dimensionless. We assume that 
at this stage our Lagrangian possesses a CP-symmetry and that these couplings are 
all real. All CP-violation will come from spontaneous symmetry breaking. The fields 
which get VEVs are 0i, 02, 02, A ^3R, and u c . The first 4 break D 3 symmetry, 
while (A^r) = A® R T3R breaks SU(2)k symmetry. All of these VEVs we take to be 
real. u c and u c are fields which get VEVs in the right-handed neutrino direction. 
These VEVs we take to be complex, and furthermore take the combination to be in 
a particular direction: {u) c uf) oc ce ldl + X R e ld2 . More will be said later about this 
choice of direction. All of the above mentioned VEV values are taken to be of order M* 
rather than M c , but are allowed to lie enough below M* to be able to give the desired 
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hierarchies in the Yukawa matrices. We also assume that the fields listed here, which 
get VEVs, obtain mass near the cutoff scale M*. This lifts these fields high enough that 
they cannot adversely affect the gauge coupling unification results we wish to preserve. 

Our superpotential allows 9 extra U(l) symmetries which we take to be symmetries 
of the theory in order to forbid problematic superpotential terms. 6 We choose to 
parametrize these symmetries by the following 9 fields: 1P3, ip^, ^a, A A 3 r, 

u c , and u c . We choose each of these fields to have an arbitrary charge under one 
and only one U(l). After specifying these charges, the charges of the remaining fields 
are determined by the terms in the superpotential. In addition to these U(l)s, the 
superpotential can have Z2 symmetries which are extensions of the first orbifold Zi- 
The transformations of the bulk fields under this symmetry have already been given. To 
have such a symmetry, we can choose A 3 r odd and the remaining 7 independent fields 
(H+_ is already determined) even under this Z%. The transformations of the remaining 
fields can then be determined from the superpotential terms. Another possible choice 
is ipa odd with the remaining fields even. It can be shown that it is not possible to 
similarly extend the second Z2 orbifold to the brane terms of the superpotential. 

Our 5D bulk Higgs field H(x^, y) can be represented as a 4D massless state and a 
KK tower of 4D massive states. In order to find the effective Yukawa matrices, we need 
to know the overlap between the 5D bulk field H(x^,y) and the effective 4D massless 
Higgs H°(x^). 7 This overlap is 



H(x^ y) D e^nnp^fMlH ^) (3.13) 



where 



"v = \/^TT (3 - 14) 



9 



2M C 



The kink mass dependence gives the massless Higgs field an exponential profile, local- 
izing it to an end of the extra dimension. We will use n x , n H , n L , n R , etc. to stand 
for a normalization as above with the particular value for the X quantum number sub- 
stituted for the X in equation ( 3.14|) . In the same manner, at times we will use n\ or 



6 We can avoid unwanted flavor changing neutral currents from the breaking of these extra symme- 
tries by instead assuming only Z n subgroups of the U(l)s sufficient to forbid the unwanted superpo- 
tential terms. 

7 When addressing the overlap between 5D and 4D fields, to avoid confusion we will explicitly show 
the dependence of the fields on the spacetime coordinates. 
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n_ 3 to refer to these normalization factors with the indicated X value inserted. We 
can see explicitly in equation (|3.13| ) that the up- and down-Higgs wavef unctions are 



localized towards opposite ends of the y direction domain [0, ^1 due to their opposite 
X quantum numbers. 

The light fields which will correspond to our observed three families of particles 
come from various mixings of the fields in the superpotential. Consider the following 
terms contained in W%: 

W 2 D pV*. (a^a + bo^K)] 8{y - ^f) (3.15) 
With the factors of M* explicit: 

&„(^) V ^(c e * + ^) = & (3.16) 

One combination of ip 1 ^ and N£ is massive while the other combination leads to a 
massless field. In terms of a and b, the overlap between the massless field ipa° an d the 
original fields is: 

€ 3 n c r a ° (3-17) 
N c a D -~n c ^° 

1 



1 _L _2i 
1 + |6|a 



The light fields ip° form a D3 doublet and correspond to the 1st and 2nd family right- 

b 



handed particles. The factor of | in the overlap between N% and will lead to 



l ^ ce' 9 i+x e w 2 * n ^ e (2,2) elements of the Yukawa matrices. This factor is important 
in fitting inter-family fermion mass ratios and in providing us with nontrivial phases in 
the Yukawa matrices. 

Now we turn to a different subset of W<i- 

W 2 D TE a {d y - m x )16 a (3.18) 
+ [2<7^16 a 16 3 ] 5{y) + [2 V (V++) a Xa] 5{y - — ) 

These mass terms lead to mixing between the brane fields 163(0;^), Xa{%^) and the bulk 
fields lQ a (x tl ,y). The overlap between the resulting left-handed massless field (ip^(x f1 )) 
and those fields in the superpotential: 

Mx") 3 ^ 3 V) (3-19) 
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Xa{x' ) D — e 2 



The overlap in the right-handed fields is 



(3.20) 



The definitions of and £ (y) have been given in equation fl3.7|). ^3 and ■0 : 

correspond to the left- and right-handed 3rd family fields. Other definitions follow: 

1 



cO 



n R = 



1 + r 2 
1 



-4 + 4e x ^ 



(3.21) 



1 + r 



2 1 



cr 



r = 



P 1 



+ 



For more information on our treatment of these overlaps and normalizations, please see 
Appendix |B[ 

The left-handed 1st and 2nd family fields are equal to ip a , which do not mix with 
any other fields. We replace all fields within W\ with their massless components. 
This yields the X-dependent Yukawa matrices for the massless fields with left-handed 
doublets on the left: 



2T, 



( "oiS e j 



Y 



\ 



V 



'A 
b 





- an HL Bn^— -Yn^ 
"° LR ^° b LR "0 L 



A. 



(3.22) 



Definitions follow:* 



L = e 2 ni, 

R = e 2 

A = \zpn L n R n H 



(3.23) 



3 We use the fact that for all Yukawa terms, Xl + Xr + Xh = 0. 
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_ K (A) 
A 



A) 

7o 



A 3 M* b' 

02 

£o = 7o — 

01 

The Yukawa matrices may be simplified by looking at the normalization constant 
rf. We assume that a <C \b\, and so we can approximate n c ~ 1. This is not incompat- 
ible with our definition of /3 , as we also expect (3 <C 1. 

An approximation may also be made within n^. 

n L = , 1 (3.24) 



I +r 2 



1 +r 2 ^ 

"I 

We assume that i] > (9(1), and since p ~ 1/40 for reasonable values of the compact- 
ification and cutoff scales derived from gauge coupling unification, ^e XhC ' <C ^- for 
— 10 < ( < 10, which easily encompasses the C-range which has a chance of fitting the 
data. 

It is convenient to reparametrize r 2 in terms of other variables. With the definition 
r 2 may be rewritten as 

r 2 = j 2 {ll + el). (3.26) 
With this redefinition the normalization constants ni and hr become 

nx = , ', . , (3.27) 



i + 4^ 

P 2 n 2 x 



It is this definition for the nx which we will use throughout the rest of this paper. 
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Under the approximation of rf ~ 1, the Yukawa matrices simplify to: 



/ 



Y 



i i ? r r 

n l a b 'o l 2r 3fl 
77f? Po-T"T7? 70 -F- 



V 



o 



6 Li? 




\ 



/ 



A. 



(3.28) 



These are the Yukawa matrices which we will analyze in the next section. 9 

The lowest order diagrams which contribute to the Yukawa matrices follow. For 
each diagram we give the Yukawa element (s) to which it contributes. 

H 

(3,3) 



V>3 



H 



A 



3R 



(t>a 



;i,3) (2,3) 



H 



(2,2) 



H 



A 



(1,2) (2,1) 



4. Analysis 

We take two routes in the analysis of our model. In the first subsection we use ana- 
lytic methods to extract predictions from our theory. Relations between the Yukawa 
parameters and the observables are given but are not solved due to the complexity of 
the equations. A more precise numerical analysis, in which a full fit is achieved, is then 
given in the following subsection. 

9 We could have included terms in our superpotential allowing (3,1) and (3,2) elements in the 
Yukawa matrices. However, as long as such matrix elements are hierarchical, they cannot affect the 
theoretical values of the low energy observables. Such terms are also not necessary in our theory (e.g., 
we do not have a left-right symmetry which requires their presence) and so we have left them out of 
our superpotential. 
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4.1 Analytic Fitting 



The starting point for our analysis is the Yukawa matrix of equation (|3.28 ). In our 
model, this Yukawa matrix is defined at the compactification scale M c . We have taken 
into account the effects of integrating out the heavy fields but have neglected the 
running effects on the Yukawa matrix elements between M* (where the superpotential 
of equations ( |3.11| ) and Q3.12| ) is defined) and M c . We assume that these running effects 
are small enough to ignore at the level of a few percent. 

To make the analysis easier, we define a u , a^, etc. below. 



Y(X) 



/ 







■«n^ — 

u n L n R 



ao^r^ — 
u n L n R 



A) 



ce iy l +X R e t0 2 n L n R 








After inputting the proper X quantum numbers: 



3R- 



7o2T ; 



3R- 



-X L C/2 
~X L C/2 



X x 



(4.1) 



a, 



(m) 2 



1 



lu = -7oe 
e u = — Eqc 



-C/2J_ 

ni 

-C/2J_ 
ni 



a d = ane^^ — 

Pd = ce iH -3e lB 2 P° e ^ m^iL 

ld = 7oe" 
e d = e e' 

X b 



-C/2J_ 

ni 

-C/2J_ 

ni 



X t e-^ 



= a e^^4 — 



Pe = 

7e = 7oe 
£ e = e e 
A 



i 



ce w l +e i9 2 

= ^p3C/2^1_ 
"-3 

= :! C/2^ 

n_3 



A>e c ~1r 



(4.2) 



and 



1 



1 + 



/■■ J 7g+eg 



(4.3) 



There are eleven parameters associated with the fermion masses derived from the 
Yukawa matrices: (, a , /3q, 70, £o ; c > 81, @2, K , X t , and tan/3. Because there are 13 
independent observables in the quark and charged lepton sectors: 9 masses, 3 quark 
mixing angles, and the CP-violating phase in the CKM matrix, we can expect two 
predictions in this model. 

We assume that both of the quark Yukawa matrices are hierarchical ffl]. As we 
explained within |^6| this allows us to use a simple set of rotations to diagonalize our 
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quark mass matrices and leads to a simple CKM matrix: 

- Sl2 - S f 3S * 3 1 S* 23 + S V 2 S* 13 J. (1.1) 

S 12 S 23 — s 13 ~ s 23 — s 12* s 13 1 

where 



(4.5) 



The eigenvalues of the diagonalized quark Yukawa matrices lead to the quark masses: 

- A *V2 mt — m c - |/3„| 2 



*12 


_ Pu ' 


6 12 


~ a d 
_ p d > 


Sl2 


— *12 


6 125 


„J7 
b \Z 


~ rr 
C M! 


5 13 




■513 


— 6 13 


6 13' 


JJ 
b 23 


~ 'V 
/U) 


6 23 


- 7d, 


■S23 


— 5 23 


- S U 
b 23- 



^2 m 6 — m a — |/3 d | 



(4.6) 



We will take the charged lepton Yukawa matrix to be lopsided to some extent, with 
7 e ~ 0(1) and the rest of the matrix following a hierarchy. 10 This choice leads to the 
charged lepton masses: 

We want to use the kink mass, here parametrized by C, to localize the two Higgs wave 
functions to opposite branes. As stated before, if H u is localized towards the SO(10) 
brane and is localized towards the PS brane, then there is a natural reason why m t 
is larger than m& while m u is smaller than m^. As can be seen from the a x (equation 
( |4.2| )), which play a large part in determining the masses of the first family fermions, 
if ( > we have m u suppressed by and enhanced by e^. Similarly, the third 
family masses have the opposite dependence as shown in the relation: = \ t e~^^^-. 
Looking at j x , we see that ( > enhances 7 e while suppressing ^ u and 7^, so it is not 
unreasonable to assume that 7 e ~ 1 while 7«,7d <C 1. 

With a few approximations within the CKM matrix the mixing angles and CP- 
violating angle (3 are: 



\Vus\ - \S* 12 + S 13* S 23| - |Sl2 



did 



_ a u (3 d 



(4- 



V cb\ ^ 1^3 + -Sr2 s 13l - I ^23 1 ^ \jd ~ 7«l 



10 



The lopsided effect with order one 7 e is crucial to achieve correct b-r unification 
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VI I * U * * I I I 

i ub\ — | S 13 S 12 S 23\ — 



1 



13 = arg -tt-^ ~ arg 1 - - arg 1 

With hierarchical structure in the up and down quark matrices it is the case that 



W, 



cb\ 



(2,3) p _ (2,3)u 
(3,3) D (3,3)u 



Here this is reflected in the relation \V c b\ — |7d — 7«|- To have 
a nonzero V c b, we need to have some dependence on right-handed quantum numbers 
within 7 X . We have found that in our setup, regardless of the placement of fields, the 
kink mass effects (which come with (X) can only bring X L into the Yukawa element 
ratio For this reason we have added the field A 3R to the superpotential, and it 

is the VEV of this field (proportional to T 3R ) which gives us the necessary difference 
between 7^ and 7 M to generate a nonzero V c b- 

The structure of the (2,2) Yukawa elements, proportional to ce ^ 1 ei e 2 , helps us 
to fit several observables. The Xr dependence gives different phases between the up 
and down elements, both of which contribute to the CKM matrix CP-violating phase (3. 
The phase of (3 U is also important in fitting the size of \V u b\. In addition, the magnitude 
of these (2, 2) elements helps us to fit the first to second and second to third family 
mass ratios. 

We will now show the predictions of this model in the form of relations between 
observables at the compactification scale. Consider the fermion mass combinations 
and (mA/fmA. 



\m c ) V"V/ 



— ) (M c ) ~ ^a/iT^ = v/TT^f (4.9) 



One prediction of this model is then 



\m c J \m^J \m T J 
Consider also the determinants of the down quark and charged lepton mass matrices. 

m d m s m b (M c ) = | det Y d \ (^j = a 2 d X b (^j (4.11) 



mem^m^Mc) = \ detY e \ ( -j=j = a e \ T y-j= 



\ 3 / \ 3 

Vd\ 2^ I V d 
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Because a d = a e and \ b = A T , we have the (exact) prediction: 



m d m s m b (M c ) = m e m^m T (M c ) (4.12) 

These two predictions hold regardless of the actual values taken on by the Yukawa 
parameters. 

As for the Yukawa parameter values themselves, these must be determined from 
the eleven remaining independent observables, if possible. The dependence of these 
observables on the parameters is complicated enough that we find it impossible to 
make a fit nonnumerically. Actual fit values for the Yukawa parameters must wait 
until the numerical analysis in Section [4.2| . 

Our predictions are relations at the compactification scale, while our data has been 
taken at the weak scale. In order to determine the extent to which our predictions 
are reasonable, we need to estimate the renormalization effects on the masses between 
the two scales. In this analysis, we choose to diagonalize the Yukawa matrices at the 



compactification scale, and use the RG formalism of Barger, Berger, and Ohmann [58 
to relate the observables at M z (or at m t in the case of the top quark) to their values 
at M c with simple scaling relations. 



M diag (M c 



'S u (M z ,M c )m u (M z ) 

S u (M z ,M c )m c (M z ) J ( 1.13) 

S t {m u M c )m t {m t ) 



M d diag (M c 



S d (M z ,M c )m d (M z ) 

S d (M z ,M c )m s (M z ) 

S b (M Zl M c )m b (M z ) 

M e diag (M c ) = 

-S e (M z ,M c )m e (M z ) 

S e (M z ,M c ) mfM (M z ) 

S T (M z ,M c )m T (M z ) 

|Vf (M c ) = 

\V ud \ 2 {M z ) \V US \ 2 {M Z ) S(M z ,M c )\V ub \ 2 (M z ] 

\V cd \ 2 (M z ) \V CS \ 2 {M Z ) S(M z ,M c )\V cb \ 2 (M z ) 

S{M z ,M c )\V td \ 2 {M z ) S(M z ,M c )\V ts \ 2 (M z ) \V tb \ 2 {M z ) 
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The scaling factors S can be found in [ |3| . In deriving these scaling factors, the authors 
have used 2-loop running and have included running effects from the gauge couplings 
and the third family Yukawa couplings. We make a further approximation in that we 
keep only 1-loop effects and neglect all running effects from the Yukawa sector except 
for the top Yukawa coupling. With these approximations 

(4.14) 



S U {M Z . 


M c ) 


~y~%m t ,M c )G u (M z ,M c ) 


S t (m t 


M c ) 


~ yf(m t ,M c )G u (m t) M c ) 


S d (M z 


M c ) 


~G d (M z ,M e ) 


S b (M z 


M c ) 


^y^(m t ,M c )G d (M z ,M c ) 


S e (M z 


M c ) 


~ G e (M z ,M c ) 


S T (M Z 


M c ) 


~G e (M Z} M c ) 


S(M Z 


M c ) 


~ y 2 t {m t ,M c ) 



with 



y t {m u M c ) = exp 
G X (M Z , M c ) = exp 



, , A t 2 0)d(log/x) 



167T 2 

a,(M c ) 



mt 
M c 

M z 



(4.15) 



^^(^d^Og/i) 



(4.16) 



The c x and 5 govern the MSSM 1-loop running effects due to the gauge couplings: 

(4.17) 



H 3 , — 

15' ' 3 



B 



7 16 

15' d ' 3 



9 

?, 3, 



1, -3 



We estimate that the neglect of the running due to and A T introduces 5% error 
in the down and charged lepton sectors and around 1% error in the up sector. We 
further estimate the neglect of 2-loop gauge running to be ~ 1% error in G u and Gd, 
and ~ 4% error in G e . 
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Because we use MSSM running from M c down to Mz, errors are introduced when 
we include particles in the running below their mass scales. These effects can be 
taken into account by weak-scale threshold corrections. However, in our analysis we do 
not keep track of the supersymmetric particle masses and so these corrections cannot 
be calculated. We can, nevertheless, estimate these effects. The SUSY thresholds 



for all of the down sector quarks were calculated in | 50| . There are tan /3-enhanced 
diagrams which contribute to m b , m s and m^. Although the Higgsino contribution to 
the d, s threshold corrections are Yukawa-suppressed, the dominant correction comes 
from gluino and Wino loop diagrams which affect all three of these quarks equally 
(up to some small differences due to unequal squark masses). We therefore make the 
simplifying approximation that the SUSY threshold corrections for all of the down 
quarks are the same. The same approximation can be made separately for the up and 
charged lepton sectors where the differences in thresholds among the same particle type 
arise only from squark mass differences, and are small. 

By assuming some range of SUSY parameters, Pierce et. al. [[19| have estimated the 
SUSY threshold corrections to the 3rd family masses. Assuming a value of tan/3 ~ 30 
and \i > 0, we estimate from plots in ^{| the SUSY threshold effects (at m t for top and 



Mz for bottom and tau) in terms of % shifts and errors. The shifts for these particles, 
and by approximation the shifts for the 1st and 2nd family particles, are 

m u ,m c ,m t (+2.0 ±3)% F t = 1.02 

m d ,m s ,m b (-10.0 ± 10)% F b = 0.90 (4.18) 
m e ,TO M ,m T (+1.5 ±2)% F T = 1.015. 

The factors F x will be used to keep track of these corrections. 

SUSY threshold corrections to the CKM matrix elements were also calculated in 
50| . As those authors state, to a good approximation the only CKM elements which 



shift their sizes are \V ub \, \V cb \, \V t d\, and \V ts \, and these shifts are the opposite of the 
chargino-induced shift in m b . In addition, J shifts approximately twice the amount 
of the CKM elements. The contribution from charginos to the mj shift is plotted 



separately in Pierce et. al. [49|, and we estimate it to be 5 ± 10%. For those CKM 

(4.19) 



observables which have large SUSY threshold corrections: 

\V ub \, \V cb \, \V td \, \V ts \ (-5.0 ± 10)% F v = 0.95 



J (-10.0 + 20)% Fj ee 0.90. 



As explained in | 5"0|| , each side of the unitarity triangle has one element which has a 
large SUSY threshold. This leads to the threshold for J, and also implies that the 
angles of the triangle are unaffected. Thus sin 2(3 does not get a correction. We assume 
that the correction to \ek\ is approximately the same size as the correction for J. 
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The mutual dependence of yt(m t , M c ) and X t means that there is no simple analytic 
way of integrating to find yt{m t , M c ). We choose, therefore, to fit m t (m t ) by adjusting 
X t (M c ) and using the 1-loop equations with X t and gi only to run down to m t . The 
1-loop RG equation is: 



~dt 



Xt 
16tt 2 



(4.20) 



with t = log /i where \i is the energy scale. The relation between the experimental value 
of mAmA and XAmt): 



v v 
m t (m t ) = X t (m t )F t —sm(3 ~ X t {m t )F t -j= 



(4.21) 



We have already assumed tan/? ~ 30, which means that sin/? is negligible, and we have 
included the thresholds for m t (m t ) in the scale factor F t . 

The PS-breaking VEVs of x° an d X c are °f order the cutoff scale. As in f2"D|| , we 
choose to parameterize these VEVs with a dimensionless parameter Cbrane^ 



(x c ) = (r) 



AM* 



^"fi f 5 Cbrane 



(4.22) 



Naive dimensional analysis leads to a value for Cbrane : 

Cbrane 0.27 



(4.23) 



Using this value and using the assumption of 5D gauge unification as in |20j with the 
4D GUT scale inputs 



Mq = (2.5 ±0.5) x 10 16 GeV 
« GUT = 1/(24 ±1) 
e 3 = -0.035 ± 0.005 



(4.24) 



leads to knowledge of M*, M c , and ai(M c ): 





= 2.3 x 10 17 GeV 


M c 


= 2.4 x 10 14 GeV 


«i(M c ) 


= 0.035 


a 2 (M c ) 


= 0.040 


a 3 (M c ) 


= 0.044 



(4.25) 
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The uncertainties in the 4D GUT scale parameters (especially q;gut) introduce ±9% 
errors into G u and G d and ±2% error into G e and y t . The errors listed here on G u and 
Gd have a high correlation and can be mostly neglected when considering ratios. 

Considering all the sources of errors discussed above, we choose to parametrize 
our theoretical errors as ±9% G u and Gd (with the errors highly correlated), ±4.5% 
on G e , and ±2% on y t . In addition, in our expressions, the masses themselves should 
have extra theoretical errors from both threshold corrections and the neglect of the 
contribution to Yukawa running due to A& and A T . We choose the up-type quarks to 
have ±3% error, down-type ±11% error, and the charged leptons ±5%. Each of these 
errors is correlated to some degree within each particle type. 

Starting from the «j(M c ) and using 1-loop running, we find the gauge coupling 
scale factors to be 

G U (M Z , M c ) = 0.33 ± 0.03 (4.26) 
G d (M z , M c ) = 0.33 ± 0.03 
G e (M z ,M c ) = 0.70 ±0.03. 

Fitting the central value of m t (m t ) = 169 ± 4GeV, allowing for the m t la range, and 
including the 3% threshold error on the top quark leads to 

A t (M c ) = 0.59 ± 0.12 (4.27) 
y t (m t ,M c ) = 0.903 ±0.018. 

We have included the 2% theoretical error on y t . 

We are now in a position to analyze our predictions. First, consider 

(^v(^W)-(-W < 4 - 28 > 

\m c J \"W \m T J 

We choose to make a prediction for (jf^j (M z )- The prediction and corresponding 
experimental values are: 



\ m cJ th 



69 

\ " l c / cxp 



( ( mb Fb Gd 1 

\m^J \m T F T G e y t 
0.0037 ±0.0006 



(M z ) (4.29) 



0.0023 ± 0.0010 

cxp 



Our prediction falls outside of the experimental bounds and leads to a 1.2cr discrepancy 
with the data. Our model favors a larger value for j , indicating possibly a larger 
m u and/or smaller m c than currently measured. 



- 28 - 



Consider the second prediction: 



m d m s m b (M c ) = m e m^m T (M c ) (4.30) 
After using the scale factors to get down to M z , this becomes 



{m d m s m b )th{M z ) 



3 



(M z ) (4.31) 



~ (10.7 ± 5.0) x 10 5 MeV 3 
(m d m s m b ) exp (M z ) = (5.2 ± 2.6) x 10 5 MeV 3 . 

The uncertainty of 50% in the predicted value comes mostly from the 9% errors in the 
gauge running scale factors G e and G d and the 11% error associated with the down-type 
masses. The experimental error is dominated by the uncertainty in m s (and hence m d ). 
The discrepancy in the two values is about la, and we predict a slightly larger scale 
for the down-type masses than measured. 

4.2 Numerical Fitting 

Within this section we apply numerical methods to fit our model to the data. By using 
automated techniques, we can find all of the Yukawa parameters which fit the data, 
including those parameters which were difficult to determine analytically. We can also 
easily include 2-loop gauge running and (1-loop) running due to the whole Yukawa 
matrices, not just X t . 

Our fitting procedure starts from the Yukawa matrix in equation ( |3.28| ) taken at 
the compactification scale. We determine the gauge couplings at that scale assuming 
unification and some usual values for the unification parameters as listed in equations 
). Using 2-loop gauge and 1-loop Yukawa MSSM running, the Yukawa matrices are 



run from the compactification scale down to Mz and diagonalized to find the fermion 
masses and quark mixing angles. On the way to Mz, the top running mass m t (mt) 
and the running due to Xt below m t (yt(Mz,m t )) are determined. The observables at 
M z are shifted by SUSY threshold effects, listed in equations ( [4.18|) and (|4.19|) . To 



counteract the inclusion of the incorrect top running between Mz and m t , we multiply 
the observables at Mz by appropriate powers of y t (Mz,m t ). The correct powers are 
determined from the scale factors listed in the previous section equations ( |4.14| ). The 
observables are then compared to data in a x 2 function, which we minimize by altering 
the Yukawa input parameters. 

The sources of error in this section are the following. First, as in the analytic 
section, the threshold effects on the observables at the weak scale due to the SUSY 
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Table 4: Observables used in the x 2 analysis. The theoretical errors are combinations of 
estimates of weak threshold effects, GUT parameter uncertainties, and in the case of ek an 
additional theoretical uncertainty from Bk- All observables are at the Mz energy scale, 
except for the top mass, which is at mt- 
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±UlcLl H/11U1 /0 
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9^ 1 


o.u 


9^ A 


n 70 _i_ n i 7 HpV 


/vyi 1 /vyi 1 

rrityrrit) 
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9 


a n 


lUy 31 1 VjC V 


rn 3 
m d 


8 5 

O . 'J 


2 5 


8 8 


18 9 + 1 7 


rn s 
m b 


24.8 


2.2 


24.9 


0.0199 ±0.0049 


m b 


2.1 


12.2 


12.4 


2.91 ±0.36 GeV 


m e 
«l M 


0.0003 


2.0 


2.0 


0.004738 ± 0.000095 


m M 
m T 


0.02 


2.0 


2.0 


0.0588 ±0.0012 


m T 


0.02 


2.0 


2.0 


1.747 ±0.035 GeV 


\v us \ 


1.6 


0.0 


1.6 


0.2240 ±0.0036 


\Vcb\ 


1.9 


10.0 


10.2 


0.0415 ± 0.0042 


\v ub /v cb \ 


9.3 


0.0 


9.3 


0.086 ±0.008 


\Vtd\ 


9.8 


10.0 


14.0 


0.0082 ±0.0011 


sin 2(3 


6.5 


0.0 


6.5 


0.739 ± 0.048 


J x 10 5 


10.0 


20.0 


22.4 


3.0 ±0.7 


\ £ k\ 


0.7 


26.6 


26.6 


0.00228 ±0.00061 



particle spectrum are not calculated. We estimate these effects in the same way as 
in the previous section in equations ( |4.18| ) and ( |4.19| ). Second, the uncertainties in 
our canonical 4D GUT parameters still introduce theoretical errors in our low energy 
data. By altering these inputs and studying their effects on the calculated observables 
numerically, we estimate the theoretical errors on our observables from these effects. 
These estimates are in general a bit less than the errors assigned in the analytic sec- 
tion because some correlations are automatically included numerically which were not 
included before. Third, we have the usual experimental errors. We choose to combine 
errors in quadrature and assign a combined theoretical and experimental error to each 
observable. These observables and errors are listed in Table |j. Most theoretical errors 
cancel in the same-family mass ratios, up to a few percent. Hence, we have assigned a 
minimum 2% error on these ratios due to threshold effects. This is especially important 
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Table 5: Observables, target values, best fit values, and x 2 contributions. All observables 
are at the Mz energy scale, except for the top mass, which is at wi[. 



Observable 


Target Value 


Fit Value 


X 2 Contribution 


Q 


22.7 ± 1.7 


23.5 


0.21 


m c 


0.73 ± 0.17 GeV 


0.59 GeV 


0.67 


m t (m t ) 


169 ± 7 GeV 


167 GeV 


0.11 


m a 
m d 
m 3 
m b 


18.9 ± 1.7 


17.0 


1.28 


0.0199 ±0.0049 


0.0238 


0.64 


m b 


2.91 ±0.36 GeV 


2.46 GeV 


1.57 


m e 

TTln 


0.004738 ± 0.000095 


0.004729 


0.01 


m M 
m T 


0.0588 ±0.0012 


0.0588 


0.00 


m T 


1.747 ±0.035 GeV 


1.757 GeV 


0.08 


\v us \ 


0.2240 ± 0.0036 


0.2237 


0.01 


\v cb \ 


0.0415 ± 0.0042 


0.0412 


0.00 


\v ub /v cb \ 


0.086 ±0.008 


0.090 


0.21 


\Vtd\ 


0.0082 ±0.0011 


0.0084 


0.03 


sin 2/5 


0.739 ± 0.048 


0.720 


0.15 


J x 10 5 


3.0 ±0.7 


3.0 


0.00 


\ £ k\ 


0.00228 ±0.00061 


0.00210 


0.09 






Total: 


5.06 



for the lepton ratios, whose experimental errors are negligible. 

Taken from Martin and Vaughn |]5"1"| |, we use 2- loop gauge and 1-loop Yukawa 



renormalization group equations. As explained before, our boundary conditions are 
at M c : Given a set of choices of the Yukawa parameters, we have the 3x3 complex 
Yukawa matrices at M c ; with the assumption of gauge unification and with the 4D GUT 
parameters (equation fl4.24|) ), we have aj(M c ). Running due to the neutrino sector of 
the theory is neglected. We fit the 11 Yukawa parameters to the 16 observables listed 
in Table (| by minimizing a \ 2 function 

X 2 = J2 Ilv ■ ( 4 - 32 ) 



X i 



Our fit to the data is in Table [5] and the corresponding Yukawa parameters for 
this fit can be found in Table [(]. A true y 2 function assumes gaussian errors for its 
observables, something which we have implicitly assumed but which is not true for 
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Table 6: Minimum x 2 fit Yukawa parameters. 



Parameter 


Value 


c 


2.152 


Q'o 


0.0002007 


A) 


0.003274 


7o 


0.02187 




0.0009816 


c 


1.492 


Oi 


5.090 




2.718 


K 


1.495 


A t (M c ) 


0.6057 


tan /3 


25.42 



some of the observables used. Our x 2 is more of an indication of how good the fit is, 
and the minimization of the function is a method by which we can find a "best" set of 
parameters for the fit. The fit value for our x 2 is around 5, with the majority of the 
contribution coming from the down-type quark sector. 

The value for m u /m c in the numerical fit (0.0040) is about 0.5<r away from the value 
found in the analytic section (0.0037 ± 0.0006), and is consistent with our prediction 
of a larger m u and/or smaller m c than measured. In addition, the best fit value for 
m^m^mfe in the x 2 analysis (5.0 x 10 5 MeV 3 ) is consistent with the value from the 
previous section ((10.7 ± 5.2) x 10 5 MeV 3 ) (within 1.1a) and with the measured value 
((5.2 ± 2.6) x 10 5 MeV). As for the free parameters, in the analytic section we fit only 
one of these: X t (M c ) = 0.59±0.12. The numerical fit value found here, 0.6057, falls well 
within the range of the analytic fit value. Thus, the results of the numerical (Section 
OP and analytic (Section |4.1|) analyses are consistent. 



5. Justification in 6D 

In this section we discuss a 6D version of our theory which naturally justifies our setup 
given in this paper. Our 5D analysis has assumed that the PS brane keeps its symmetry 
even though x c and x c get their VEVs near the cutoff scale and break the PS symmetry 
entirely. As mentioned in section ||], there are generic large corrections in the Kahler 
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potentials of PS-brane fields ifj 

This problem can be solved in a geometric way by the addition of a sixth dimension 
along which \ c an d t/> are separated. Even for a tiny sixth dimension, I/R2 ~ M*/3, PS 
breaking effects are suppressed by e~ 3 which is enough suppression to keep our Yukawa 
relations. All of the numerical analysis given in the paper remains the same as long as 
the sixth dimension is not too large compared to the cutoff scale. 

There are two distinct ways of constructing 6D models. First, a 6D N=l (4D N=2) 
theory with T 2 /Z2 gives four fixed points (3+1 dimensional spacetime). All the fields 
in the bulk in the 5D theory are now in 6D in the 6D theory, while all the fields living 
on branes still live on branes. As we have one 10 dimensional hypermultiplet and two 
16 dimensional hypermultiplets, we need one additional 16 dimensional hypermultiplet 
in order to cancel the 6D irreducible gauge anomaly of SO(10) f52fl . 

Second, a 6D N=l theory with T 2 /(Z2 x Z' 2 ) also gives a 4D N=l theory below 
the compactification scale. There are four fixed lines along x 5 = 0, x 5 = ttRi/2, x 6 = 0, 
and x 6 = ttR 2 /2 and four fixed points at the corners. The gauge sector is extended to 
6D, but all of the other bulk fields can still be on 5D fixed lines. This setup is anomaly 
free without introducing additional states. 

Gauge coupling unification restricts the possible sizes of R\ and R2 in both cases. 
However, we can easily recover the 5D theory used in this paper in the limit I/R2 — > M*. 
Note, as long as I/R2 is close to, but somewhat less than, M* we simultaneously have 
gauge coupling unification (as described in the 5D formulation) and PS symmetry 
relations for Yukawa couplings. 

We choose to focus on a 6D N=l theory with T 2 /Z2. Such a theory has been 
studied by Asaka, Buchmuller and Covi [jl9[ where the SO(10) gauge group in the 



bulk is broken down to SM x U(l) x by two different Wilson lines. One breaks SO(10) 
down to the Pati-Salam gauge group along x 5 , and the other breaks SO(10) down to 
SU(5) x U(l) x along xq. Therefore, we end up with four fixed points; SO(10) brane at 
{x^^xq) = (0,0), Pati-Salam brane at ( 2L f 1 ,0), Georgi-Glashow brane (SU(5) x U(l) x ) 
at (0, ^) and flipped SU(5) brane (SU(5)' x U(l)^) at sfa). We have provided a 
schematic of the extra-dimensional space in Figure |l|. The 5D setup we have considered 
so far can be easily lifted up to this configuration. The procedure is the following. 

• 5D bulk states are extended to 6D bulk states. Extra hypermultiplets are intro- 
duced to cancel the 6D anomaly. 

• Fields on the 5D SO(10) brane are located on the 6D SO(10) brane. 
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SU(5)xU(l)x 


SU(5)'xU(l)' x 


SO(10) 

i 


PS 



e [0, ?f] 

Figure 1: A representation of the 6D orbifold space. The orbifold fixed points are at the 
corners of the space. At each corner a different subgroup of SO(10) is preserved. X5 runs 
along the horizontal direction from to ^y 1 , while xq runs along the vertical direction from 
to 2f^. 



• Fields on the 5D Pati-Salam brane are located on the 6D Pati-Salam brane except 
X c and x c - Note, in this framework we no longer need x° an d X'\ since the 
gauge group is broken down to 811(3)^ x SU(2) L x U(l) y x U(l) x by the two 
Wilson lines. Nevertheless, any PS breaking effects on the SO(10) or PS branes, 
emanating from the PS breaking fixed points, are suppressed. U(l)x is broken 
near the compactification scale by 16kmk and X c kink5 the fields which generate the 
nonzero kink mass. 

The 6D model is simple and economical. The additional states introduced to cancel 
the 6D anomaly are hypermultiplets and can become heavy by themselves. In addition 
they do not affect the differential running of the three gauge couplings. Let us focus on 
gauge coupling unification. The spectrum of massive Kaluza-Klein vector multiplets 
are given by 

M + - = v © 2 + ^) 2 
M- + =A/(^) 2 + (f) 2 



(SMxU(l) x ) 
PS 



SM x U(l) x 
/SU(5)\ 
V SM J 
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SO(10) 



Figure 2: A representation of the SO(10) group space. The space has been broken into four 
subspaces by the orbifold symmetry, under which each subspace has a different parity. 



M__ 



// 2n+l N 2 I / 2m+l \2 ( SU(5) X U(l) x 

V 1 *i ) +l R * ] V SMxU(l) x 



The breakdown of SO(10) into these subgroups is illustrated in Figure |2]. If we send 
i?2 — > 0, we recover a 5D theory in which SU(5) x U(l) x is broken to the SM x U(l) x . 
The R2 — > limit gives the same result as in Hall and Nomura (and also in Kim 
and Raby gj). Thus we fix M c ~ l/R x at around 10 14 GeV with M* ~ 10 17 GeV. If 



i?2^* ~ 2 or 3, we get a tiny correction from extra — h states and the result would 
be proportional to |(6ps — &sm) l°g(^") which is neglegibly small. 11 The elongated 6D 
rectangular configuration gives a perfect setup for the construction of a realistic SO(10) 
model. 

The only remaining question is the applicability of bulk field localization in 6D 
which was possible with a kink mass in 5D. Though the quantitative results of 6D 
localization are different from the 5D case, all of the qualitative aspects remain the 
same 



53|. More interestingly, 4D N=l supersymmetry is preserved. 



L The notation used is defined in Kim and Raby |2C 
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6. Summary and Discussion 

Extra dimensions provide a nice framework for understanding how grand unified theo- 
ries may be realized in nature. In this paper we have constructed a 5D SO(10) model 
which accommodates the quark and charged lepton masses and the CKM matrix. The 
model uses 11 parameters to fit the 13 independent observables of the quark and charged 
lepton Yukawa sectors, allowing us to make two predictions: ^(Mx) = 0.0037±0.0006 
and m d m s m b (M z ) = (10.7 ±5.0) x 10 5 MeV 3 , both of which are roughly la larger than 
the experimental values. The kink mass localizes the bulk Higgs fields according to 
their U(l)x quantum numbers, giving some explanation for the hierarchy m& <C m t 
and m d > m u . Our 5D SO(10) model can be considered to be an effective theory com- 
ing from 6D SO(10) with one small and one large extra dimension. If the size of the 
sixth dimension is very small (i.e. if the inverse of its characteristic length is near the 
cutoff scale) then all of the calculations done here for the 5D model can be regarded as 
good approximations for the 6D case. 

In a theory in which there are additional (heavy) vector-like states, commonly called 
Froggatt-Nielsen (FN) fields, which have the same quantum numbers as the light states, 
the light and heavy states can mix. The extra dimension naturally gives us a chance 
to unify Froggatt-Nielsen fields with ordinary matter fields if both come from the same 
hypermultiplet in higher dimensions. In our model, however, we chose to place the 1st 
and 3rd family fields on opposite branes to easily take advantage of the effects of the 
kink mass, and so our model does not achieve the unification of the massless and FN 
states into bulk fields. Nevertheless, it may be possible to place all of the matter fields 
in the bulk and still restrict the Yukawa terms of those fields to opposite branes to gain 
the desired hierarchy from the kink mass. We leave this possibility for further research. 

Further work is needed in six main areas in order to make our model complete. 
First, the neutrino sector should be included to explain neutrino oscillation exper- 
iments. We have already fixed the neutrino Dirac masses in terms of the Yukawa 
matrix between the left and right-handed neutrinos, but the heavy Majorana masses of 
the right-handed neutrinos have not yet been determined. It would be very interesting 
to expand the model to include Majorana masses and to investigate the resulting neu- 
trino masses and mixings. Second, the weak scale supersymmetry breaking mechanism 
should be specified. As there is as yet no such universally accepted mechanism, we 
chose not to specify how SUSY is broken in our model. As a consequence, we could not 
calculate the electroweak threshold corrections resulting from superparticle spectrum. 
Extra dimensions provide new interesting channels for the understanding of supersym- 
metry breaking and its mediation and can give entirely different superparticle spectra. 
We leave the weak scale supersymmetry breaking physics to future work. Third, our 
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analysis is somewhat incomplete since we have taken values for ogut, and £3 from 
other sources. Hence our treatment of gauge and Yukawa coupling RG running is not 
completely self consistent. We have however accounted for this shortcoming by includ- 
ing a theoretical uncertainty obtained by varying the gauge coupling parameters at Mq- 
While a more complete unification treatment would be preferable, this would require 
knowledge of supersymmetry breaking and the sparticle spectrum. Fourth, we have 
neglected the effects of running on the Yukawa matrices between the compactification 
and cutoff scales. Our assumption is that such effects would contribute at the level 
of a few percent. Further research is necessary to investigate these effects. Fifth, the 
UV completion of the higher dimensional gauge theory could give us a better under- 
standing of the model. String theory does not allow arbitrary matter configurations 
and the constraints coming from it are usually stronger than those from field theory. 
Therefore, it would be interesting to see if the 5D model considered here can be derived 
as an effective field theory from a string theoretic starting point [54. 55]. Finally, we 
have assumed that many fields in our theory obtain vacuum expectation values at cer- 
tain scales and in particular directions in group space. A complete theory would need 
to justify these VEVs by the minimization of the appropriate potentials. Such work 
would be most easily done after finding a successful UV completion for the theory, as 
many of these VEVs are of the order the cutoff energy scale. 
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Appendix 

A. D3 Family Symmetry 

In this section, we define the D 3 group, give its character table, and give other infor- 
mation necessary to understand the couplings and representations used in our model. 
Our presentation is an abbreviation of a more complete description given in . 



D 3 is the group of all rotations in three dimensions which leave an equilateral 
triangle invariant. The group contains six elements in three classes: E; C3, Cf; C a , 
Cf,, C c . E is the identity element. C3 and Cf signify 120° and 240° rotations about an 
axis through the center of and perpendicular to the triangle. C a , Cb, and C c signify 
180° rotations about the three different axes which run from the center to the three 
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Figure 3: Graphic representation of our choice of D3 group elements. The solid (black) lines 
represent the triangle. The 3 dashed (blue) lines represent the axes about which we take the 
C a , Cb, and C c rotations. The center (red) circle and dot represent the axis about which we 
take the C3 rotation (it points out of the page). Our rotation convention is right-handed. 



vertices of the triangle. We choose orientations such that Cb = C a Cs and C c = C a C\. 
Figure [| contains a graphic representation of our conventions for the group elements. 
The group has two inequivalent singlet representations: 1a and 1b, and one doublet 
representation: 2a- 

The character table is 



D 3 


E 


c 3 


c a 


1a 


1 


1 


1 


1b 


1 


1 


-1 


2a 


2 


-1 






(A.1) 



We choose our representations as follows. When acting on a one dimensional rep- 
resentation, the elements are the characters in the character table. When acting on the 
two dimensional representation, the elements are 



E 



1 
1 



e 
e- 1 



C a 



1 

1 



(A.2) 



with 



2m/3 



. The remaining elements can be found by group multiplication. 



Our Lagrangian of Section I!2" contains terms with various combinations of D3 
fields. It is understood that the Lagrangian only contains the D 3 singlet (1a) part of 
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these combinations. Here we list the combinations of fields we have used, and their 
singlet projections. 

Let Qi be 2 a (doublet) fields under D3, and let the internal degrees of freedom be 

represented by Xi and y,i\ 0j = ( X% ) . Also, let be a 1b ("anti-symmetric" singlet) 

\ViJ 

field, with internal degree of freedom a: (f) = a. Then the combinations we need are: 



Q\ x 6*2 1 i A = X\Vi + IJ1X2 (A. 3) 

0i x 2 x 0| 1a = {xxy-2 - yix 2 )a 
0i x 2 x 3 | = x ± x 2 x 3 + y±y 2 y 3 



Of course, multiplication by a "true" singlet 1a has no effect concerning the D 3 structure 
of the Lagrangian terms. 

B. Massless States and Wavefunctions 

This section will show how we determine the overlap between the original Lagrangian 
states and the massless states. We choose to illustrate this through a specific example: 
the determination of the overlap between the Lagrangian fields and the 3rd family 
left-handed states in equations (|3.19|) of section 

The relevant superpotential terms are (from equation ( |3.12| )): 



W 2 D ^ ++ ) a (d y -m x )(^) a (B.l) 

— — wE 

+ [2a^ a %) + [2i7 (^ ++ ) a Xa\ S(y - — ) 

The equations of motion are: 
dW 

7u7-^ = (B ' 2) 

= (d y - m x ) (V-) a + 25(y)a^-y^M'^ 3 + 25(y - ^)r /v ^ Xa . 

We have added factors of where necessary to ensure that we have unitless couplings 
a and 77. This equation is satisfied for the massless projections of these fields. Away 

from the branes, (ip ) oc e mxy . Because (ip ) is odd under both orbifold parities 

and because e mxV is not an odd function about either y — or y — the profile 

of (ip ) a must have discontinuities at both branes. These discontinuities can cancel 

the brane terms in the equations of motion. Define the overlap between ip 3 and the 
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massless field ipt as ^3 ^ n L^t- We can then solve for the other two fields: 



(B.3) 



Xa D -e 2 



What remains is to determine the normalization n^. We do this by requiring that the 
effective 4D field ip® have a canonical Kahler potential term. Listing the fields which 
contribute to this term (we have suppressed the gauge field factors here): 



ttR 

n L a a J 



(B.4) 



The integral is 



2 

d?/ e ^ 



-R 

■2 



2X L Cy a 

ay e * R 



(B.5) 



which implies that for canonical normalization: 



ELM 



n L r] z 



Choose n L to be real: 



n L 



1 + r 



4 + 4e x ^ 

n L rf 



(B.6) 



(B.7) 



The other mixings and normalizations listed in Section |3.2| are calculated in the 
same way. 



C. A Lifted 4D Model 

We present here a 5D version of a 4D model by Dermisek and Raby in [32] which itself 
was based on prior works |56|, f)7], [58], j59[. Our purpose is to illustrate how such a 4D 
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model could be placed into a 5D context and what advantages can be gained from the 
use of the extra direction. 

In using the extra dimension, we have separated the fields on the PS brane from 
the important mass parameter M x = Mo(l + aX) which gives much of the distinction 
between the particle types. There are matter fields on the PS brane, and bulk matter 
fields which mediate between these brane fields and the SO (10) M x VEV in such a way 
to give us the desired Yukawa matrix elements. 

The setup for this model is largely the same as that presented in Section |3|, save 
for the following points. This model does not have a kink mass in the bulk; the small 
up quark mass comes instead from an approximate left-right symmetry in the model. 
Here we use a B — L VEV on the PS brane instead of T^r. Other differences lie in the 
placement of the matter fields and the extra VEV fields. As with the model presented 
in the main body of the paper, for the current model under discussion we also have in 
mind a setup which includes a small 6th dimension (as detailed in Section ||) in order to 
separate fields on the PS-brane from the PS-breaking VEVs. This allows us to preserve 
the PS Yukawa relations on the PS-brane below the cutoff scale. 

The states are placed as follows. In the bulk, we have 8 matter hypermultiplets, 
forming 4 doublets under D3 and transforming as 16s under SO(10). Each of these 
hypermultiplet doublets has a different parity under the orbifold. The Higgs, as before, 
is contained inside a 10 hypermultiplet of SO(10). These fields are listed in Table [7]. 



On the SO(10) brane there is only the mass parameter M x , which is a singlet under D3 
and is a mix of singlet and X under SO(10): M x = X x (*^) = M (l + aX). On the 



Pati-Salam brane, we have three sets of left- and right-handed matter fields. Two of 
these sets form a doublet under D 3 . In addition, there are several extra fields: a , a , 
A, Axs, $1,, and Of these extra fields, those with subscript a are D 3 doublets. The 
rest are 1a save A which is 1b under D 3 . All of these extra fields get nonzero VEVs 
except for <pi. All extra fields are SO(10) singlets, save Ai§ which is an SU(4) adjoint 
and gets a VEV (A l5 ) = (A° l5 ) (B - L). The PS brane fields are listed in Table g. 
We choose the following superpotential: 



Wt = % - —) {X^H^l + \^ a H + A 2 ^ ++ ) a H^ L ) (C.l) 

W 2 = W a dyl6 a + W a dyl6' a + lE a d y 16 a + I6' a dyl6' a + 10 dy 1 




nR 




+6(y 




+- 



+ 



) 



a L 



\ 3 A 15 cj) a ilj 3 + \ 4 A 15 (j) a il) a + \ b Ai}) a $ L 
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Table 7: 



Bulk fields 



Field PS Symm D 3 Symm 




There are 12 U(l) symmetries associated with our superpotential. We choose to 
parametrize these by allowing the following 12 fields to be charged each under different 

U(l)s: Vs, H c + _, M x , (V-) a , (^-+) a > (^-+) a ' $L ' Al5 ' After s P ecif y in g 
the U(l) symmetries of the above fields, the U(l) transformations of all other fields 
are uniquely defined. There are in addition 2 Z 2 symmetries. First is the Z 2 orbifold 
parity under y — > —y. The transformations of the bulk fields under this symmetry have 
already been defined. In addition, let the rest of the independent fields: ip 3 , ^3, M x , 
$ L , A 15 be even under this symmetry. The second Z 2 involves a sign ambiguity in the 
transformation of Under a given symmetry, if $l — > e w $L, then the superpotential 
terms imply that $_r — > e m7r e ta &R with n — 0, 1. We choose to require that n — 1 here 
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Table 8: PS Brane fields 



Field PS Symm D 3 Symm 





(4,2,1) 


2a 




(4,1,2) 


2a 


^3 


(4,2,1) 


1a 




(4,1,2) 


1a 





(1,1,1) 


2a 





(1,1,1) 


2 a 


A 


(1,1,1) 


1b 


A 15 


(1,1,3) 


1a 




(1,1,1) 


1a 




(1,1,1) 


1a 



in order to forbid terms created by the replacement of $l by or vice- versa. We also 
assume that /M* C 1 so that replacements like <£>z, — > $fj are negligible. 12 Let 
the 12 independent fields be uncharged under this symmetry. We require all of these 
symmetries just listed to be symmetries of the theory so as to forbid unwanted extra 
terms in the superpotential. 

The superpotential has a left-right symmetry, under which tp 3 <-> ip^, ip <-> ip c , 

if) „ (E\ (B 

16 J \1&J \16 
with the family D 3 symmetry, is broken spontaneously by the VEVs of and 
which we require to be slightly different. This difference is encapsulated by a small 
parameter 77: ($_r) = (3>l) (l + 77). 

We now turn to the equations of motion for the left-handed states: 

dW 

= 0^d y ^ ++ ) a = (C.2) 




This symmetry, which commutes 



,y x = o a, M ) =0 

S = 0„ + 5(y)M x ty/ + _) a = 



12 Such terms given by the substitution $l — > <&fj or $^ — > $| would lead to a Yukawa matrix 
structure different from the one desired and so should be forbidden by some symmetry. 
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dW 



d y U> 



d y U'__)+5{y)M x ^ ++ ), 



S(y - ^) A ;i .l.r/. t + A i A i -/ < ; I ' , -■ A-, .4 I '2 



$r = 



dW 



d y U'_)+5(y)M x ^ ++ ) 1 



+S(y 



i 

nR 



^3^1501^3 + A 4 A 15 2 V ; 2 - A 5 AV>1 



$r =o 



Solving these equations leads to knowledge of the overlap between the massless fields 
and the original fields in the Lagrangian. The following relationships only have the 
massless components on the right hand sides of the equations. (We have replaced the 
brane fields by their VEVs) 



^) [A3(ALi 5 )(0i)^3 + A 4 (A1 15 )(0 2 )^ 2 -A 5 (A)^i <<I>l) 
(ft-) 2 D ^(y) [A 3 (A 15 > (0 2 )^ 3 + A 5 {A) fa] ($ L ) 



(C.3) 



(^++)i 



M, 



A 3 (Ais) ^3 + A 4 (A 15 ) ^0 2 ) ^ 2 - A 5 (A) VlJ ($L 

[A 3 (A 15 > (0 2 )^3 + A 5 {A) fa] ($ L ) 



The equations for the right-handed states can be obtained from the left-right sym- 
metry present in the model. We list the most important of these equations: 



1 



\ 3 (A 15 )( ( f )l )r 3 + \ 4 (A 15 )U 2 )r2-^(A)^ (C.4) 



« + ) 2 [^3 (A 1B > (0 2 > ^ + As (A) n (*r) 



The three ipi fields span the space of the left-handed massless states. Because the 
other fields which we've integrated out have massless components, the kinetic energy 
and gauge interaction terms for the ipi fields are no longer orthonormal. The same is 
true for the ijjf states by the left-right symmetry. We will discuss the effects of rotating 
and rescaling these fields to an orthonormal basis later. 
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We replace the (ip ++ ) a and \ip++) fields in W\ by their corresponding massless 
parts in order to get the low energy Yukawa matrices. The result: 13 

/ e'p renT u c \ 
Y u = | -e'p ep reT u c ] A (C.5) 

K = 



Definitions follow for these variables, where we have used M x = M (l + aX), 
(®r) — (1 + v)i an d have added factors of the cutoff scale in order to make the 
couplings Xi all unitless. We have also assumed C 1 and a ~ 0(1). Tf represents 
the (B — L) quantum number for the field /. 

A 2 A 4 (W0 2 W> 2 4a 

£ = \ / _ cq 6) 

3AiM3M (l + a)(l-3a) V ' ; 

A 2 A 5 (A) ($ L > 2 4a 




e' 



XiM^M (l + a)(l-3a) 
277(1 - 3a) 



4a 

A 3 (0i) 3(1 -3a) 



a = 



A 4 

1 + a 
1 — 3a 



4a 



v !2M r 



These Yukawa matrices are the same as those found in | 32|j , except for the (1,3) and 



(3,1) elements. It has been shown in [46] that (1,3) elements are needed in models of 
this kind in order to fit sin 2(3. 



13 We have chosen to use the notation found in [p2[ to ease comparison between prior works and our 
own. 
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Here we list the lowest-order diagrams which give the Yukawa matrices. Each 
diagram is followed by the element (s) to which it contributes. 
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H 



V>3 



(3,3) 



H 



3£- 



.4 



15 



(1,3) (2,3) 



.4 



15 



'03 



M x 
-X 



H 



(3,1) (3,2) 



H 



4>a 



X- 



.4 



0a 



15 



ft 



(2,2) 



15 



M x 



(2,2) 



M x 
5C- 



.4 



ft 



(1,2) (2,1) 



.4 



V'o 



M x 



(V>++) £ 



(1,2) (2,1) 



We performed a x 2 analysis, the same as that used in [SB], save that here we have 
nonzero (1,3) and (3,1) terms and a new parameter k. The fit parameters follow. 
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First the GUT parameters: 

1 



«GUT 

la 

e 3 = -3.61% 



25.12 (C.7) 



Mn = 2.54 x 10 16 GeV 



Next the (modulus) Yukawa sector: 



A = 0.70 (C.8) 

r = 24.8 

a = 1.19 

£ = 0.0090 

p = 0.061 

e' = 0.0034 

k = 0.15 

Next the Yukawa phase information (cf) x = arg(x)) in radians: 

(f) a = 0.51 (C.9) 
P = -1.87 
<J) K = 0.89 

The other Yukawa parameters are assumed to be real. Finally some SUSY breaking 
scales and tan f3: 



H = 234 GeV (C.10) 
Mi /2 = 606 GeV 
m 16 = 4160 GeV 



A = -7736 GeV 

2 

= 1.80 



= 0.128 
tan/3 = 52.2 



The fit itself can be found in Table |^. As explained in Section T?2 ) we are using 
our x 2 function as a vehicle to find the best possible fit rather than as a true statistical 
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Table 9: Best fit. Mass dimensions are in GeV. 



Observable 


Target Value 


Fit Value 


X 2 Contribution 


1 


137.04 ±0.14 


137.0 


0.08 


C*EM 








x 10 5 


1.1664 ±0.0012 


1.166 


0.11 


oc s 


0.1172 ±0.0020 


0.1167 


0.06 


M Top 


178.0 ±4.3 


176.5 


0.12 


m b 


4.220 ±0.090 


4.243 


0.07 


M b — M c 


3.40 ±0.20 


3.35 


0.06 


m s 


0.089 ±0.011 


0.104 


1.86 


^xlO 3 


2.03 ±0.20 


2.00 


0.02 


m d 
m s 


0.050 ±0.015 


0.074 


2.54 


M T 


1.7770 ±0.0018 


1.777 


0.00 




0.10566 ±0.00011 


0.1057 


0.13 


M e x 10 3 


0.51100 ±0.00051 


0.5110 


0.00 


V us 


0.2230 ±0.0040 


0.2216 


0.12 


v cb 


0.0402 ±0.0019 


0.0390 


0.40 


Vub 

v cb 


0.0860 ±0.0080 


0.0863 


0.00 


e K x 10 3 


2.28 ±0.23 


2.35 


0.10 


M z 


91.188 ±0.091 


91.19 


0.00 


M w 


80.419 ±0.080 


80.41 


0.01 


m c 


1.30 ±0.15 


1.15 


1.00 


(b -> 57) x 10 3 


0.334 ±0.038 


0.335 


0.00 


sin 2f3 


0.727 ±0.036 


0.700 


0.57 


V td x 10 3 


8.20 ±0.82 


8.35 


0.03 






Total: 


7.30 



X 2 function. Our best x 2 is 7.30, indicating that we are not fitting some of the data. 
Specifically, as shown in Table |9|, the observables m s , and m c have x 2 contributions 
of 1 or greater and make up the majority of the x 2 fit value. The fit values of the 
down and strange quark masses are on the large side, while the charm mass fit value 
is smaller than the data. We present this fit as a first attempt at this kind of model. 
More work is needed to alter the model to obtain a better fit. 

In our analysis we have used a basis in which the massless matter fields are not 
orthonormal. Rotation and rescaling to a canonical orthonormal basis would in general 
introduce changes to the Yukawa matrices. We have neglected effects from this change 
of basis, and our justification for this is the following. Were a fit to be done with these 
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effects included, the input Yukawa parameters would compensate by changing their 
values. We assume that the input parameters could compensate to the extent that we 
would obtain essentially the same fit in this case as in the case we have presented here 
in the appendix without these extra effects. We leave it to further research to explore 
whether this assumption is a good one. 

References 

[1] S. Dimopoulos, S. Raby and F. Wilczek, Phys. Rev. D 24, 1681 (1981). 

[2] S. Dimopoulos and H. Georgi, Nucl. Phys. B 193, 150 (1981). 

[3] L. E. Ibanez and G. G. Ross, Phys. Lett. B 105, 439 (1981). 

[4] N. Sakai, Z. Phys. C 11, 153 (1981). 

[5] M. B. Einhorn and D. R. T. Jones, Nucl. Phys. B 196, 475 (1982). 

[6] W. J. Marciano and G. Senjanovic, Phys. Rev. D 25, 3092 (1982). 

[7] S. Raby, Phys. Rev. D 66, 010001 (2002). 

[8] H. Georgi and C. Jarlskog, Phys. Lett. B 86, 297 (1979). 

[9] Y. Kawamura, Prog. Theor. Phys. 105, 999 (2001) [arXiv:hep-ph/0012125]. 

[10] Y. Kawamura, Prog. Theor. Phys. 105, 691 (2001) [arXiv:hep-ph/0012352]. 

[11] G. Altarelli and F. Feruglio, Phys. Lett. B 511, 257 (2001) [arXiv:hep-ph/0102301]. 

[12] L. J. Hall and Y. Nomura, Phys. Rev. D 64, 055003 (2001) [arXiv:hep-ph/0103125]. 

[13] L. J. Hall and Y. Nomura, Phys. Rev. D 65, 125012 (2002) [arXiv:hep-ph/0111068]. 

[14] L. J. Hall and Y. Nomura, Phys. Rev. D 66, 075004 (2002) [arXiv:hep-ph/0205067]. 

[15] A. B. Kobakhidze, Phys. Lett. B 514, 131 (2001) [arXiv:hep-ph/0102323]. 

[16] A. Hebecker and J. March-Russell, Nucl. Phys. B 613, 3 (2001) 
[arXiv:hep-ph/0106166]. 

[17] A. Hebecker and J. March-Russell, Nucl. Phys. B 625, 128 (2002) 
[arXiv:hep-ph/0107039]. 

[18] R. Dermisek and A. Mafi, Phys. Rev. D 65, 055002 (2002) [arXiv:hep-ph/0108139]. 

[19] T. Asaka, W. Buchmuller and L. Covi, Phys. Lett. B 523, 199 (2001) 
[arXiv:hep-ph/0108021]. 



- 50 - 



[20] H. D. Kim and S. Raby, JHEP 0301, 056 (2003) [arXiv:hep-ph/0212348]. 
[21] S. B. Giddings and A. Stromingcr, Nucl. Phys. B 307, 854 (1988). 
[22] S. R. Coleman, Nucl. Phys. B 310, 643 (1988). 
[23] G. Gilbert, Nucl. Phys. B 328, 159 (1989). 

[24] Y. Kawamura, H. Murayama and M. Yamaguchi, Phys. Rev. D 51, 1337 (1995) 
[arXiv:hep-ph /9406245] . 

[25] K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 83, 2522 (1999) 
[arXiv:hep-ph/9906271]. 

[26] L. M. Krauss and F. Wilczek, Phys. Rev. Lett. 62, 1221 (1989). 

[27] C. H. Albright, Int. J. Mod. Phys. A 18, 3947 (2003) [arXiv:hep-ph/02 12090]. 

[28] M. C. Chen and K. T. Mahanthappa, Int. J. Mod. Phys. A 18, 5819 (2003) 
[arXiv:hep-ph/0305088]. 

[29] G. Altarelli and F. Feruglio, arXiv:hep-ph/0405048. 

[30] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D 61, 033005 (2000) 
[arXiv:hep-ph/9903417]. 

[31] H. D. Kim, J. E. Kim and H. M. Lee, Eur. Phys. J. C 24, 159 (2002) 
[arXiv:hep-ph/0112094]. 

[32] R. Dermisek and S. Raby, Phys. Rev. D 62, 015007 (2000) [arXiv:hep-ph/9911275]. 

[33] M. Battaglia et al, arXiv:hep-ph/0304132. 

[34] A. Hocker, H. Lacker, S. Laplace and F. Le Diberder, Eur. Phys. J. C 21, 225 (2001) 
[arXiv:hep-ph/0104062]. For updated Lepton-Photon 2003 result, see: 
http: //www. slac.stanford.edu/ laplace/ckmfitter/ckm_results_summer03.html. 

[35] K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66, 010001 
(2002). 

[36] H. Leutwyler, Phys. Lett. B 378, 313 (1996) [arXiv:hep-ph/9602366]. 

[37] R. Gupta and K. Maltman, Int. J. Mod. Phys. A 16S1B, 591 (2001) 
[arXiv:hep-ph/0101132]. 

[38] C. W. Bauer, Z. Ligeti, M. Luke and A. V. Manohar, Phys. Rev. D 67, 054012 (2003) 
[arXiv:hep-ph/0210027]. 

- 51 - 



[39] P. Azzi et al. [CDF Collaboration], arXiv:hep-ex/0404010. 

[40] H. Fusaoka and Y. Koide, Phys. Rev. D 57, 3986 (1998) [arXiv:hep-ph/9712201]. 

[41] T. Inami and C. S. Lim, Prog. Theor. Phys. 65, 297 (1981) [Erratum-ibid. 65, 1772 
(1981)]. 

[42] N. Arkani-Hamed, T. Gregoire and J. Wacker, JHEP 0203, 055 (2002) 
[arXiv:hep-th/0101233]. 

[43] R. Kitano and T. j. Li, Phys. Rev. D 67, 116004 (2003) [arXiv:hep-ph/0302073]. 

[44] H. Murayama and A. Pierce, Phys. Rev. D 65, 055009 (2002) [arXiv:hep-ph/0108104]. 

[45] L. J. Hall and A. Rasin, Phys. Lett. B 315, 164 (1993) [arXiv:hep-ph/9303303]. 

[46] H. D. Kim, S. Raby and L. Schradin, Phys. Rev. D 69, 092002 (2004) 
[arXiv:hep-ph/0401169]. 

[47] S. M. Barr and I. Dorsner, Phys. Lett. B 556, 185 (2003) [arXiv:hep-ph/02 11346]. 

[48] V. D. Barger, M. S. Berger and P. Ohmann, Phys. Rev. D 47, 2038 (1993) 
[arXiv:hep-ph/9210260]. 

[49] D. M. Pierce, J. A. Bagger, K. T. Matchev and R. j. Zhang, Nucl. Phys. B 491, 3 
(1997) [arXiv:hep-ph/9606211]. 

[50] T. Blazek, S. Raby and S. Pokorski, Phys. Rev. D 52, 4151 (1995) 
[arXiv:hep-ph /9504364] . 

[51] S. P. Martin and M. T. Vaughn, Phys. Rev. D 50, 2282 (1994) [arXiv:hep-ph/9311340]. 

[52] T. Asaka, W. Buchmuller and L. Covi, Nucl. Phys. B 648, 231 (2003) 
[arXiv:hep-ph/0209144]. 

[53] H. M. Lee, H. P. Nilles and M. Zucker, Nucl. Phys. B 680, 177 (2004) 
[arXiv:hep-th/0309195]. 

[54] T. Kobayashi, S. Raby and R. J. Zhang, arXiv:hep-ph/0403065. 

[55] T. Kobayashi, S. Raby and R. J. Zhang, "Searching for realistic 4d string models with 
a Pati-Salam symmetry: Orbifold arXiv:hep-ph/0409098. 

[56] T. Blazek, S. Raby and K. Tobe, Phys. Rev. D 60, 113001 (1999) 
[arXiv:hep-ph/9903340]. 

[57] R. Barbieri, L. J. Hall and A. Romanino, Phys. Lett. B 401, 47 (1997) 
[arXiv:hep-ph/9702315]. 



- 52 - 



[58] R. Barbieri, L. J. Hall, S. Raby and A. Romanino, Nucl. Phys. B 493, 3 (1997) 
[arXiv:hep-ph/9610449]. 

[59] T. Blazek, M. Carena, S. Raby and C. E. M. Wagner, Phys. Rev. D 56, 6919 (1997) 
[arXiv:hep-ph/9611217]. 



- 53 - 



